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Abstract 



We study long-time dynamics of a class of abstract second order in time evolution equations 
in a Hilbert space with the damping term depending both on displacement and velocity. This 
damping represents the nonlinear strong dissipation phenomenon perturbed with relatively 
compact terms. Our main result states the existence of a compact finite dimensional attractor. 
We study properties of this attractor. We also establish the existence of a fractal exponential 
attractor and give the conditions that guarantee the existence of a finite number of determining 
functionals. In the case when the set of equilibria is finite and hyperbolic we show that every 
trajectory is attracted by some equilibrium with exponential rate. Our arguments involve 
a recently developed method based on the "compensated" compactness and quasi-stability 
estimates. As an application we consider the nonlinear Kirchhoff, Karman and Berger plate 
models with different types of boundary conditions and strong damping terms. Our results 
can be also applied to the nonlinear wave equations. 

AMS 2010 subject classification: Primary 37L30; Secondary 37L15, 74K20, 74B20 
Keywords: Nonlinear plate models; state-dependent damping; global attractor; dimension. 

1 Introduction 

We study a class of plate models with the strong nonlinear damping which abstract form is the 
following Cauchy problem in a separable Hilbert space H: 



We impose the following set of hypotheses: 

Assumption 1.1 (A) The operator A is a linear self-adjoint positive operator densely defined 
on a separable Hilbert space H (we denote by | • | and (•, •) the norm and the scalar product 
in this space). We assume that the resolvent of A is compact in H. We also denote by H s 
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(with s > 0) the domain V(A S ^ 2 ) equipped with the graph norm | • \ s = \A s ^ 2 -\. In this case 
H- s denotes the completion of H with respect to the norm | • \_ s = \A~ S / 2 • |. Below we 
denote by {e^} the orthonormal basis in H consisting of the eigenfunctions of the operator 
A : 

Aet = Afcefe, < Ai < A 2 < • • • , lim A^ = oo, 

k— >oo 

and by Pjv the orthoprojector onto Spanje^ : = 1,2,..., N}. 

(D) For some value 9 G (0, 1] the damping operator D maps H\ x Hq into H_g and possesses 
the properties: 

(i) For every g > there exist a Q > and f3 s > such that 

(D(u,v),v) > a g \v\j and \D(u, v)\_ e < P e \v\e (2) 

for all (u;v) G Hi x Hq, \u\ 2 + \v\ 2 < g 2 (it is allowed that a e — > and (3 e — > oo as 
g — >■ oo). 

(ii) For every £> > there exist 7 e > and C e > such that 

(ZV V) - I>(u 2 V) V - v 2 ) > 7 ,|^ - v 2 \ 2 - C e \u l ~ u 2 \U{l + \v l \ 2 e + | V 2 |i) (3) 

for all (u l ;v l ) £ Hi x Hq, \u % \\ + \v l \ 2 < g 2 , and for some 5 > (it is allowed that 
7 e — > and C e — >■ oo as £ — > oo). 
(hi) We assume that for every u G Loo(0,T;i7i) the mapping w i— >■ D(u,v) is weakly 
continuous from L2(0,T; Hq) into £2(0, T; H^q) and 

|Z>(u\«) - £>(u 2 >«)l-i < C.lu 1 - u 2 \i- s {l + Mi) 
for all I it* 1 1 , 1 17 1 < g, u l £ Hi, v G -ffg, and for some 5 > 0. 
(iv) In the case 9 < 1 we have 

^(u 1 V) - D{u 2 ,v 2 )U < C e [\v l - v 2 \ e + In 1 - u 2 |i(l + \v x \ e + \v\)] (4) 

for all (u l ;v l ) G Hi x iZg, + \v l \ 2 < g 2 . In the case 9 = 1 we assume a stronger 
(compared to (4) with 6 = 1) inequality: 

\Diu\v 1 ) - D(u 2 ,v 2 )\-i < C e [ft; 1 - v\ + \u x - u\_ & (\ + \v\ + |v 2 |i)] (5) 

for all (u l ;v % ) G Hi x ju l | 2 + |u l | 2 < p 2 , where (5 > 0. 

We note that the conditions in (2) follows from (3) and (4) or (5) provided D(u, 0) = for 
every u G Hi. 

(F) There exists 5 > such that the nonlinear operator F maps His into H_q and is locally 
Lipschitz, i.e., 

\F{ui) - F{u 2 )\_ e < L{g)\ui-u 2 \i-5, ^\A l,2 Ul \<g, (6) 

We also assume that F(u) = U'(u), where U(u) is a (7 1 -functional on Hi = V(A 1 ^ 2 ), and ' 
stands for the Frechet derivative. We assume that II(u) is locally bounded on Hi and there 
exist 77 < 1/2 and C > such that 

v\A 1/2 u\ 2 + U(u) + C > , u e Hi = V(A 1/2 ) . (7) 
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Remark 1.2 Our conditions can be relaxed in different directions. For instance, for the well- 
posedness in Theorem 2.1 we can assume that 5 = in (3) and (6), and require F to map H\- a 
into H_i continuously for some a > and / > 0. In a similar way, instead of (Diii) we can 
assume other continuity properties of the damping operator D that will allow us to perform the 
limit transition in the corresponding Galerkin approximations. To obtain the global attractor 
existence, we can also relax (3) by including into its latter term the expression elu 1 — u 2 \\ with a 
small parameter e (see, e.g., Assumption 3.21 in [10] for a similar requirement in the case of the 
monotone damping). Moreover, instead of (3) we can assume that 



(D(u\ v 1 ) - D(u 2 ,v 2 ),v' - v 2 ) > le \J -v 2 \ 2 - - u 2 \ 2 [e + C 8 {e) flt^ + \v 2 \ 2 )] (8) 



for every e > 0. However, we do not pursue these possible generalizations because our abstract 
hypotheses are motivated by the plate models described below. We note that in this paper we 
concentrate on the case when 6 is positive. For some results for the case when 9 = and D(u, ut) 
is linear with respect to ut, we refer to [8]. It is also worth mentioning that our damping operator 
D is positive (see the first relation in (2)) but not monotone (see (3)) in general. Thus, we cannot 
apply here the theory developed in [10] . 

Our main applications are plate models (with hinged boundary conditions, for definiteness) . 
In this case 

• A = (— A/)) 2 , where Ad is the Laplace operator in a bounded smooth domain $7 in M 2 with 
the Dirichlet boundary conditions. We have then that H = Li{£l) and 



We also have that H s = (H 2s n for 1/2 < s < 1 and H s = H$ s (n) for < s < 1/2. 

Here above H a (Q) is the corresponding L2-based Sobolev space, Hq(Q) is the completion 
of C£°(0) in H a (Q). 

• The damping operator D(u, ut) may have the form 



where <7o(si), &i(si, «2, S3) and g(s\, S2) are locally Lipschitz functions of Sj G R, i = 1, 2, 3, 
such that <7o(si) > 0, a±(si, S2, S3) > and g(si,S2)s2 > 0. Also the functions <j\ and g 
satisfy some growth conditions (for a more detailed discussion of properties of the damping 
functions we refer to Section 5 below). We note that every term in (9) represents a different 
type of damping mechanisms. The first one is the so-called viscoelastic Kelvin-Voight 
damping, the second one represents the structural damping and the term g(u, ut) is the 
dynamical friction (or viscous damping). We refer to [20, Chapter 3] and to the references 
therein for a discussion of stability properties caused by each type of the damping terms in 
the case of linear systems. 

• The nonlinear feedback (elastic) force F(u) may have one of the following forms (which 
represent different plate models): 

(a) Kirchhoff model: F(u) is the Nemytskii operator 



V{A) = {ue H 4 (Q) : u = Au = on dtt} . 



D(u, u t ) = A [a (u)Au t } - div [01 (it, Vn)Vit 4 ] + g(u, u t ) 



(9) 



it 1 — y — k • div {|Vn| 9 Vit - /j,\Vu\ r Vu} + tp(u) - p(x) 



(10) 
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where k > 0, q > r > 0, u G IR are parameters, p G L 2 (P) and if G Lipi oc (IR) fulfills 
the condition 

liminf^ >-\l (11) 

\s\— >oo S 

where Ai is the first eigenvalue of the Laplacian with the Dirichlet boundary conditions. 

(b) Von Karman model: F(u) = — [u, v(u) + Fq] — p(x), where Fq G H 4 (Q) and p G £2^) 
are given functions, the von Karman bracket [u, v] is given by 

[u, v] = d 2 xi u ■ d 2 X2 v + d 2 X2 u ■ d 2 xi v - 2 • d 2 1X2 u ■ d 2 xlX2 v, 

and the Airy stress function v(u) solves the following elliptic problem 

A 2 v(u) + [u,u] = in 0, = v(u) = on dfi. (12) 

on 

Von Karman equations are well known in nonlinear elasticity and constitute a basic 
model describing nonlinear oscillations of a plate accounting for large deflections, see 
[21, 11] and the references therein. 

(c) Berger Model: In this case the feedback force has the form 



F(„) 



An — p(x), 



where k > and r € M. are parameters, p G L^VL); for some details and references see, 
e.g., [6, Chapter 4] and [10, Chapter 7]. 

Long-time dynamics of second order equations with a nonlinear damping was studied by 
many authors. We refer to [2, 15, 24, 25, 8, 18] for the case of a damping with a nonlinear 
displacement-dependent coefficient and to [9, 10, 11] and to the references therein for the velocity- 
dependent damping. Models with different types of a strong (linear) damping in wave equations 
were considered in [3, 4, 17, 23, 29], see also the literature quoted in these references. 

The main novelty of the current paper is the following: (i) we can consider the strong non- 
linear displacement- and velocity-dependent damping of a general structure (thus we cannot use 
analyticity of the corresponding model with a zero source term which takes place in the case when 
D(u,ut) = But, where B is a self-adjoint operator satisfying (2) with 6 G [1/2,1], see, e.g., [20, 
Chapter 3] and the references therein); (ii) this damping can be perturbed by low order terms. 

Our main result (see Theorem 3.1) states the existence of a compact global attractor and 
describes other asymptotic (long-term) properties of the system generated by (1). To establish 
this result we use the recently developed approach (see [9] and also [10] and [11, Chapters 7,8]). 
We first prove that the corresponding system is quasi-stable in the sense of the definition given 
in [11, Section 7.9], and then we apply general theorems on properties of quasi-stable systems. In 
the same framework we also establish a result on the rate of stabilization (see Theorem 3.2) which 
states that under some additional conditions every solution is attracted by an equilibrium with 
an exponential rate. To obtain this result we rely on some type of the observability inequality 
and use the same idea as in [10, Section 4.3] (see also [9] and [11]). 

The paper is organized as follows. In the preliminary Section 2 we discuss well-posedness 
of our abstract model and the dynamical system generation. We also recall several notions and 
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results from the theory of dissipative dynamical systems. Our main results on the global attractor 
for (1) and on asymptotic properties of individual trajectories are stated in Section 3. The proofs 
based on the quasi-stability property of the corresponding system are given in Section 4. In 
Section 5 we discuss some applications. 

Below constants denoted by the same symbol may vary from line to line. 

2 Preliminaries 

In this section we show that problem (1) generates a dynamical system. 
2.1 Well-posedness 

We first prove the existence and uniqueness of weak solutions to problem (1). We recall that a 
function u(t) is a weak solution to (1) on an interval [0, T] if 

u € Loo(0, T; V(A 1/2 )), d t u e Loo(0, T; H) n L 2 (0, T; V{A 9 ' 2 )) 

and (1) is satisfied in the sense of distributions. 

The main statement of the section is the following assertion which also contains some auxiliary 
solution properties needed for the results on asymptotic dynamics. 

Theorem 2.1 Let Assumption 1.1 be in force and (uq;ui) G % = V{A 1 ^ 2 ) x H. Then the 
following assertions hold. 

1. Problem (1) has a unique weak solution u(t) on R + . This solution belongs to the class 

W = C(R + ; V(A 1/2 )) n C\R + ; H), 
and the following energy relation 

£(u{t),d t u{t))+ !\D{u{T),d T u{T)),d T u{T))dT = £{u ,u l ) (13) 

J 

holds for every t > 0, where the energy £ is defined by the formula 

£(u , Ul ) =%,Mi) + n(u ) = 1 (M 2 + \A 1/2 u \ 2 ^j + U(u ). 
Moreover, this solution u(t) satisfies the estimate 

sup E(u(t),d t u{t))+ [ \A d/2 d t u(t)\ 2 dt <C(R) if E(u , Ul )<R 2 . (14) 
t>o J 



2. If u l {t) and u 2 (t) are two weak solutions such that E(u l (0),dtu l (Q)) < R 2 , i = 1,2, then 
their difference z(t) = n x (t) — u 2 (t) satisfies the relation 

t 

E(z(t),d t z(t)) + a R J \A^ 2 d t z(r)\ 2 dr < b R E(z(0), d t z(0))e c ^ (15) 
o 

for some constants or,6^,cr > 0. 
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Proof. To prove the existence of solutions, we use the standard Galerkin method of seeking for 
approximations of the form 

N 

u N {t) = Y J C k {t)e k , TV = 1,2,... 
k=i 

that solve the finite-dimensional projections of (1). Such solutions exist, and after multiplication 
of the corresponding projection of (1) by dtUN(t) we get that u^(t) satisfies the energy relation 
(13). By (7) we obtain that 

cqE(u ,ui) - ci < £(uq,ui) < C(R) 

whenever E{uq,u\) < R 2 . Therefore, by (2) the energy relation for uj^(t) yields estimate (14) for 
approximate solutions with the constant C(R) independent of N. Using the equation for U]\f(t) 
and also the conditions (2) and (6), it can be shown in the standard way that 

T 

\A- 1/2 d tt u N (t)\ 2 dt < C T (R), N = 1,2, ... , 
o 

for every T > 0. These a priori estimates show that (un; dtUN', 9uUn) is *-weakly compact in 

L oo (0,T;H 1 )x[L oo (0,T;H)nL 2 (0,T;H e )}xL 2 (0,T;H^ 1 ) for every T > 0. 

Thus the Aubin-Dubinsky theorem (see [26, Corollary 4]) yields that (uN',dtUN) is compact in 
(7(0, T; H\- £ x H- £ ) for every e > 0. These compactness properties make it possible to show the 
existence of weak solutions satisfying (14). For the limit transition in the nonlinear terms we use 
the property (Diii) in Assumption 1.1 and relation (6). It is also clear that 1 1-> (u(t);dtu(t)) is a 
weakly continuous function in % = V(A l l 2 ) x H. To obtain the energy relation in (13), we note 
that the function u n {t) = P n u{t) solves an equation of the form 

d tt u n + Au 11 = h(t) 

with some h G £2(0, T; H). This makes it possible to obtain a certain energy relation for u n 
which gives us (13) after the limit transition n —> 00. This also allows us to obtain the strong 
continuity properties of t ^ (u(t);dtu(t)) in H by the standard method. 
To prove (15), we note that z(t) = n x (t) — u 2 {t) solves the equation 

dttz + Diu 1 , dtu 1 ) - D(u 2 , d t u 2 ) + Az + Fiu 1 ) - F(u 2 ) = 0. (16) 

Thus, multiplying this equation by dtz and integrating from s to t we have 
t t 
E z {t) + J {D^^drU 1 ) - D(u 2 , d T u 2 ),d T z)dT = E z (s) - J (F(n x ) - F(u 2 ), d T z)dr (17) 

s s 

for any < s < t, where E z (t) = E(z(t),dtz(t)). Therefore, using (3), (6) and (14) we obtain 
that 

t t 

E z {t) + ^ j \A 9/2 d T z\ 2 dT <E Z { S )+CR j (l + ldrU^l + ldr^fylA^Z^dT, S < t, 

s s 

for some cr > 0. Now we can apply the Gronwall lemma to obtain (15) which, in particular, 
implies the uniqueness of weak solutions. □ 
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2.2 Generation of a dynamical system 

We recall that a dynamical system (see, e.g., [6, 16, 27]) is a pair (X, S(£)) of a complete metric 
space X and a family of continuous mappings S(t) t > 0, such that (i) t i-> S*(t)y 

is continuous in X for every y £ X and (ii) the semigroup property is satisfied, i.e., 5(i + r) = 
S(t) o S'(t) for any i, r > and S(0) is the identity operator. 

We also recall that the system (X, S(t)) is gradient if it possesses a strict Lyapunov function, 
i.e., there exists a continuous functional <&(y) on X such that (i) $>(S(t)y) < $(y) for all t > 
and y G X; (ii) the equality <&(y) = <&(S(t)y) may take place for all i > if only y is a stationary 
point of 5(t). 

Applying Theorem 2.1 we obtain the following assertion. 

Proposition 2.2 Let Assumption 1.1 be in force. Then problem (1) generates a dynamical sys- 
tem in the space % = V(A 1 ^ 2 ) x H with the evolution operator S{t) given by 

S(t)y = (u(t); dtu(t)), where y = (uo;ui) and u(t) solves (1). 

This system is gradient with the full energy £(uq; u±) as a strict Lyapunov function (this follows 
from the energy relation in (13)). 

We also recall that a system (X, S(t)) is called asymptotically smooth (see [16]) if for any closed 
bounded set B C X that is positively invariant (S(t)B C B) one can find a compact set K, = IC(B) 
which uniformly attracts B: sup{distx(5'(i)y, K) ■ y € B} — > as t — >■ oo. The global attractor 
(see, e.g., [1, 6, 16, 27]) of a dynamical system (X, S(t)) is defined as a bounded closed set 21 C X 
which is invariant (S*(t)2l = 21 for all t > 0) and uniformly attracts all other bounded sets: 

lim sup{distx(5'(i)y, 21) : y G B} = for any bounded set B in X. 

j>oo 

In this paper we use the following criterion of the global attractor existence for gradient systems 
(see, e.g., [28, Theorem 4.6]): 

Theorem 2.3 Let (X, S*(t)) be an asymptotically smooth gradient system such that for any 
bounded set B C X there exists r > such that J T (B) = [J t>T S(t)B is bounded. If the set 
M of stationary points is bounded, then (X, S(t)) has a compact global attractor 21 which coin- 
cides with the unstable set M. + (M) emanating from N , i.e., 21 = M + (AT). 

We recall (see, e.g., [1]) that the unstable set M. + (M) emanating from N is a subset of X such 
that for each z G M + (AA) there exists a full trajectory {y(t) : t G M.} satisfying u(0) = z and 
dist x{y{t),N) ->■ as t ->■ -oo. 

Remark 2.4 We note that we can avoid the hypothesis of Theorem 2.3 that 7 r (-B) is bounded 
if the following requirements on the corresponding Lyapunov function <E>(x) are added: (i) $(x) 
is bounded from above on any bounded set; (ii) the set &r = {x G X : <&(x) < R} is bounded for 
every R (see, e.g., [10, Corollary 2.29]). 

3 Main results 

Our first main result is the following theorem. 
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Theorem 3.1 Let Assumption 1.1 be in force. Then the dynamical system (H,S(t)) generated 
by (1) possesses a compact global attractor 21. Moreover, 

(1) 21 = M + (7V), where AT = {(u; 0) G % : Au + F(u) = 0} is the set of stationary points and 

dist H (S(t)y,Af) = inf {\S(t)y - e\ H : e G J\f} as t ->• +oo for every y^U. (18) 

(2) 77ms attractor has a finite fractal dimension. 

(3) .Any trajectory 7 = dtu(t)) :tGl} /rom i/ie attractor 21 possesses the property 

(u; d t u; d tt u) G L^R; tf 2 _ 9 x ^ x 77), (19) 

and iaere is R > such that 

sup sup (|n|!_0 + + \d a u\ 2 ) < R 2 . (20) 
7C21 teR 

(4) 77ie system (H,S(t)) possesses a (generalized) fractal exponential attractor %l exp whose di- 
mension is finite in the space Iri = Hq x H-±. 

(5) Let C = {lj : j = 1, iV} 6e a /iraie set of functionals on Hi with the completeness defect 

£c = tc{H\,H) = sup{|n| : u G Hi,lj(u) = 0, j = 1, N, \u\i < 1} . 

Then there exists eq > swc/i that under the condition eq, < £0 ^ e se ^ £ is (asymptotically) 
determining in the sense that the property 



rt+l 

lim max / \L(u l (s) — u 2 (s))\ 2 ds = for two solutions u 1 and u 2 

3 Jt 



implies that lim^oo \S(t)yi — S(t)y2\n = 0. Here above S{t)y.i = (u l (t);dtu l (t)), i = 1,2. 

We recall that the fractal dimension dimj M of a compact set M in a complete metric space X 
is defined as 

v lniV(M,e) 

dimf M = hm sup , , , , 

7 e ^o ln(l/e) 

where N(M,e) is the minimal number of closed sets of diameter 2e in X needed to cover the 
set M. 

We also recall (see, e.g., [13] and also [10, 22] and the references therein) that a compact set 
2lcxp C 77 is said to be a (generalized) fractal exponential attractor for the dynamical system 
(H,S(t)) iff 2l exp is a positively invariant set of finite fractal dimension (in some extended space 
77) and for every bounded set D C 77 there exist positive constants to, Cd and 70 such that 

d H {S{t)D |2t cxp } = supdist w (5(t)x, 2t cxp ) < C D ■ e -^- tD \ t > t D . 

As for the determining functionals, we mention that this notion goes back to the papers by Foias 
and Prodi [14] and by Ladyzhenskaya [19] for the 2D Navier-Stokes equations. For the further 
development of the theory we refer to [12] and to the survey [5] and to the references quoted 
therein (see also [6, Chap. 5]). We note that for the first time determining functionals for second 
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order (in time) evolution equations with a nonlinear damping was considered in [7], see also a 
discussion in [11, Section 8.9] We also refer to [5, 6] for a description of sets of functionals with 
a small completeness defect. Determining modes and nodes are among them. 

Using the same idea as in [9, 10, 11] we can establish the following result on convergence of 
individual solutions to equilibria with an exponential rate. 

Theorem 3.2 In addition to Assumption 1.1 we assume that F[u) is Frechet differentiable and 
its derivative F'(u) possesses the properties 

\{F f (u),w)\^ < C R Hi, w G H u (21) 

and 

\{F'(u) - F'(v),w}\-i < CrIu-v^s ■ Hi, wSH^ (22) 

for any u, v G H± such that \u\\ < R and \v\\ < R with 5 > 0. Here (F'(u),w) is the value of 
F'(u) on the element w. Let the set M be finite and all equilibria be hyperbolic in the sense that 
the equation Au + (F 1 ((f)), u) = has only a trivial solution for each (</>; 0) G TV. Then for any 
y G % there exists an equilibrium e = (</>; 0) G M and constants 7 > 0, C > such that 

\S(t)y - e\ H < Ce-V, t > 0. (23) 

We note that this type of stabilization theorems is well-known in literature for different classes of 
gradient systems, and several approaches to the question on stabilization rates are available (see, 
e.g., [1] and also [9, 10, 11] and the references therein). The approach presented in [1] relies on 
the analysis of linearized dynamics near each equilibrium and requires the hyperbolicity condition 
in a dynamical form. Here we use the method developed in [9, 10] (see also a discussion in [11]), 
and we need this condition in a weaker form. 



4 Proofs 

As it was already mentioned, the main ingredient of the proof of Theorem 3.1 is a quasi-stability 
property of the dynamical system (H,S(t)) generated by (1). 

4.1 Quasi-stability 

We show that under the conditions listed in Assumption 1.1 the system (%, S(t)) is quasi-stable in 
the sense of the definition given in [11, Section 7.9]. Namely, we prove the following proposition. 

Proposition 4.1 Let Assumption 1.1 hold. Assume that u l (t), i = 1,2 are two weak solutions 
to problem (1) with initial data yi = (uq;u\) such that \A 1 / 2 u l \ 2 + \u\\ 2 < R 2 for some R > 0. 
We denote S(t)yi = (u % (t);dtu l (t)), i = 1,2. Then there exist C(R),j(R) > such that 



\S(t) yi -S(t)y 2 \ 2 H <C{R) 



y2 \2 He -y(R)t + f e -7(K)(H| u i( T ) -u 2 (r)\ 2 dT 
Jo 



t > 0. (24) 



This type of estimates was originally introduced in [9] and related to a decomposition of the evo- 
lution operator S(t) into uniformly exponentially stable and compact parts, see also a discussion 
in [10] and [11, Section 7.9]. 

We start with two preliminary lemmas. 
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Lemma 4.2 Under Assumption 1.1 there exist T > and a constant c > independent of T 
such that for any pair u l and u 2 of weak solutions to (1) we have the following relation 



rj~i ^ rj~\ rj~\ 

TE Z (T)+ [ E z (t)dt<c\[ \z t (t)\ 2 dt+ [ \(D(t),z t 
Jo Uo Jo 



)\dt 



+ I \(D{t),z)\dt + m T (u\u 2 ) 
lo 



for every T > Tq, where z(t) = u l (t) — u 2 (t), and the functionals E z , D and are defined as 

E z (t) = E (z(t),z t (t)) = ^((z t (t),z t (t)) + (Az(t),z(t))), 
D{t) = D(u\t),u l t (t))-D(u 2 (t),u 2 (t)), 



(G(r),z t (r))dT 



+ 



(G(t),z(t))dt 



+ 



f T dt [ T (G(r),z t (r))dT 
o Jt 



with G(t) = Fiu 1 ^)) - F(u 2 (t)). 

Proof. We use the standard arguments involving the multipliers z t and z for (16). We refer 
to the proof of Lemma 3.23 in [10] and also [11, Lemma 8.3.1], where this lemma is proved 
under another set of hypotheses concerning the damping operator. However the corresponding 
argument does not depend on a structure of the damping operator. □ 

Lemma 4.3 Let u 1 and u 2 be two solutions to (1) with the initial data (v,q;u\). We assume that 

Kl 2 + Kli ^ r2 - Then 

f-T "I r f-T f-T 

max£ z (t)<c E Z (T)+ \(D(t), z t )\dt +C R / \A e,2 z t \ 2 dt+ \A l/2 z\ 2 dt , 
[°' T ] L Jo J uo Jo 

with z = u 1 — u 2 , where E z and D(t) are the same as in Lemma 1^.2. 

Proof. This follows from (17), the Lipschitz property for F and the uniform estimate (14) for 
u l (t), i = 1,2. We refer to [10, Lemma 3.25] for a similar assertion. □ 



Now we complete the proof of Proposition 4.1. Using (4) with 9 <G (0, 1] we obtain that 

\(D(t),z t )\ < C R , £ \z t \ 2 e + e\z\ 2 (l + \u\\ 2 6 + \u 2 \ 2 ) 
for any e > 0. We also have from (4) for 6 < 1 and from (5) for 6 = 1 that 

\(D(t), z)\ < C R:£ \z t \ 2 + e\z\ 2 (l + \u]\ 2 e + \u 2 \l) + CrM 2 - 
The subcritical estimate in (6) yields 



(25) 



\^t\<C r [ \z t \ 2 e dt + e ( \z\\dt + C R ^ T [ \z\ 2 dt for every e > 
Jo Jo ' Jo 
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Therefore, Lemma 4.2 implies 



TE Z (T)+ f E z {t)dt<C R , £ f \z t (t)\ 2 e dt + e f \z\\(\u\\ 2 e + \u 2 t \ 2 e )dt + C Rt ^ T f \z\ 2 dt(2Q) 
Jo Jo Jo Jo 

for every T > Tq. By Lemma 4.3 using (14) and (25) we have that 

m&xEM) < 2c E z (T) + Cr / \A e/2 z t \ 2 dt + \A 1/2 z\ 2 dt, 

[o,n Uo Jo 

where the constants cq and Cr are independent of T. Therefore, from (26) 

TE Z (T)+ [ E z (t)dt<C Rt£ [ \z t (t)\ 2 dt + e maxE z (t) [ {\u\\ 2 e + \u 2 \ 2 e )dt + C R ^ T [ \z\ 2 dt 
Jo Jo [o> T ] Jo Jo 

which, after an appropriate choice of e, implies that 

TE Z (T) + [ T E z (t)dt < C R f \z t {t)\ 2 e dt + C RtT [ T \z\ 2 dt (27) 
Jo Jo Jo 

for every T > Tq. Using (3), (6) and (17) we conclude that there exists 7r > such that 

7 R [ \z t (s)\jdt < E z (0) - E Z (T) + e [ \z{t)\ 2 dt 
Jo Jo 

+C R , e [ \z\K\ul\j + \u 2 t \ 2 e )dt + C Rt£ [ \z\ 2 dt. 
Jo Jo 

Consequently, choosing e small enough, by (27) we have that 
TE Z {T) < cr[E z (0) — E Z (T)] 

+ " 



Thus 



Cr f E z (t){\u\\ 2 + \u 2 \ 2 )dt + Cr, t [ \z\ 2 dt. 
Jo Jo 

E Z (T) < k r E z (0) + C R f E z {t){\u]\ 2 e + \u 2 \ 2 )dt + CR, T C \z\ 2 dt 

Jo Jo 



with kr < 1 and T > Tq. Now the standard argument (cf., e.g., [10, p. 62] or [11, p. 414]) leads 
to (24). This concludes the proof of Proposition 4.1. 

4.2 Completion of the proof of Theorem 3.1 

1. Proposition 4.1 means that the system (T-L,S(t)) is quasi-stable in the sense of Definition 
7.9.2 [11]. Therefore, by Proposition 7.9.4 [11] (H,S(t)) is asymptotically smooth. By Propo- 
sition 2.2 (7-L,S(t)) is a gradient system. Thus, Remark 2.4 and Theorem 2.3 imply that there 
exists a compact global attractor. By the standard results on gradient systems with compact 
attractors (see, e.g., [1, 6, 27]) we have that 21 = M + (7V) and (18) holds. 

2. Since (H,S(t)) is quasi-stable, the finiteness of the fractal dimension dim/21 follows from 
Theorem 7.9.6 [11]. 
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3. To obtain the result on regularity stated in (19) and (20), we apply Theorem 7.9.8 [11]. 

4. One can see from (1) and Theorem 2.1 that any weak solution u(t) possesses the property 

rt+l 

J {dtMr^dr <C R for all t > 0, 

provided (uo;ui) G Br = {y G H : \y\u < R}- This implies that t h-> S'(i)?/ is a l/2-H61der 
continuous function with values in % = Hg x for every y G -Br. Therefore, the existence of 
a fractal exponential attractor follows from Theorem 7.9.9 [11]. 

5. To prove the statement concerning determining functionals, we use the same idea as in the 
proof of Theorem 8.9.3 [11], see also Theorem 7.9.11 [11]. 



4.3 Proof of Theorem 3.2 

We use the same idea as in [9, 10, 11]. 

Since M is finite by (18) in Theorem 3.1 we have that for any y G H there exists an equilibrium 
e = (<fi; 0) G M such that 

\S(t)y - e\ n -> 0, t^oo. (28) 

Thus we need only to prove that S(t)y tends to e with the stated rate. 

Let S(t)y = (u(t);u t (t)). We can assume that sup t>0 \S(t)y\y, < R, for some R > 0. The 
function z(t) = u(t) — <p satisfies the following equation 

z a (t) + D(4> + z(t), z t {t)) + Az{t) + F(<f) + z(t)) - F(<f>) = 0, t > 0. (29) 

Let £{t) = E z (t) + $(£), where E z {t) is the same as in Lemma 4.2 and 

$(t) = n(0 + z(t)) - n(0) - (F(0), z) = f 1 (F(c/) + Az) - F(0), z)dA. 

■/ 

One can see that 

S(t)+ f (D(<f> + z(T),Zt(T)),Zt(T))dT = S(0). (30) 

J 

In particular, we have that £{t) is non-increasing. Moreover, since (z;zt) — > in i?i x H as 
t — )■ +00, we have that £ (£) — >■ when i — )■ +oo. Thus £(t) > for all t > 0. It is also clear from 
(6) that 

- ^(*)| < C fl |z(t)|i_i|z(t)| 9 < e|z(t)|? + C R , e |z(t)| 2 , Ve > 0. (31) 
Applying (27) for the case when u 1 ^) = u(t) and n 2 (t) = <f> we obtain that 

TE Z {T)+ [ T E z (t)dt<C R [ T \z t (t)\ 2 dt + C R , T m a x\z(t)\ 2 (32) 
jo Jo [°- T ] 

for T > To with some To > 0. Now we prove the following lemma. 
Lemma 4.4 Let z(t) be a weak solution to (29) such that 

r-T 



L 



E z (t)dt < 5 and sup E z (t) < g (33) 
T-i teR + 



12 



with some 5, g > and T > 1. Then there exists 8q > such that 



(t)\ 2 <cf \z t (t)\ 2 e alt (34) 
J o 

/or ewery < 5 < 8q, where the constant C may depend on 5, g and T . 



max \z( 

[0,T] 



Proof. Assume that (34) is not true. Then for some 6 > small enough there exists a sequence 
of solutions {z n {t)} satisfying (33) and such that 



lim < max \z(t)\ 



T n-r 

2 



n->oo I [o,T] \_Jq 

|2 



oo. (35) 



\zumdt\ 

By (33) max[ 0) T] k n (i)| 2 < f° r ah n - Thus (35) implies that 

lim / \z?(t)\jdt = 0. (36) 

n-S-oc J Q 

Therefore we can assume that there exists z* G Hi such that 

(z n ;z t n ) -»■ (z*;0) *-weakly in £«,((), T; .Hi x %). (37) 
It follows from (2) and (36) that for u n (t) = (ft + z n (t) we have the relation 

lim f \(D(u n (t),u^(t)),ift(t))\dt = for any ift G ^(0, T; 

This allows us to conclude that u* = (ft + z* G i?i solves the problem .Ait + -F(tt) = 0. From (33) 
we have that \A 1/2 {u* - <p)\ 2 < 25. If we choose 5 > such that \A 1/2 ((fti - 2 )| 2 > 26 for every 
couple (fti and <fo of stationary solutions (we can do it because the set M is finite), then we can 
conclude that u* = eft provided <5 < 6V Thus we have z* = in (37). 

Now we normalize the sequence z n by defining z n = c~ 1 z n with c n = maxr 0) yi where we 

account only for a suitable subsequence of nonzero terms in c n . It is clear from (37) with z* = 
that c n — > as n — > 00. By (35) we also have that 



T 

\z?(t)\jdt^0 as n^oo. (38) 

Relations (32) and (38) imply the following uniform estimate 

sup {|zr(t)| 2 + |^ 1/2 ^ n (t)| 2 }<C, n = l,2,.... 
te[o,T] 1 ' 

Thus we can suppose that there exists z* G Hi such that 

(z n ;z?) -»■ (z*;0) *-weakly in L^^-Hi x (39) 
The function z n satisfies the equation 

zl + —£>(</> + z n , z n t ) + Ag" + -\F(cft + z») - F(cft)} = 0. (40) 
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As above, from (2) and (38) we conclude that 

i r T 

— \(D(4> + z n ,z?),4>)\dt -> as n^oo for any tp G £2(0, T;Hg). 

c n JO 

It also follows from (21) and (22) that 

-[F(4> + z n )-F(4>)}^ (F'(<p),z*) weakly in L 2 (0,T;H^). 

Therefore, after the limit transition in (40) we conclude that z* satisfies Az* + (F'((f>),z*) = 
and, by hyperbolicity of <j) we conclude that z* = 0. Thus (39) and Aubin-Dubinski theorem [26, 
Corollary 4] imply that max^f] \z n \ — > as n — > 00, which is impossible. □ 



Completion of the proof of Theorem 3.2. By (28) we choose To > such that (33) holds 
with 5 < S and T > T . From (31) we have that £(T) < C R E Z (T). The energy relation in (30) 
and the lower bound in (2) yield that 



L 



T 

2, 



z t (t)\idt < C R £(0)-£(T) . (41) 



Therefore Lemma 4.4, relation (32) and the energy relation in (30) imply that £{T) < jr£(0) for 
some < 7h < 1. This implies that £{mT) < j%£{0) for m = 1, 2, . . . By (31), (34) and (41) we 
have that 

EjmT) < 2£(mT) + C R max \z n \ 2 < C R £(mT), m = 1, 2, . . . 

[mT,(m+l)T] 

Thus E z (mT) < Cr^r for m = 1, 2, . . . Now using (15) we obtain (23). The proof is complete. 

5 Applications 

As it was mentioned in the introduction, our main applications are plate models. 



5.1 Plate models 

For the definiteness, we concentrate on the hinged boundary conditions (the results remain true 
with other types of self-adjoint boundary conditions). Below || • || s is the norm in the Sobolev 
space H S (Q) of order s. 

Forcing term: We first check that the forcing term F satisfies Assumption 1.1(F) for all 
cases described above. 

In the case of the Kirchhoff model, the embeddings 

H%n) c L 2/( i_ a) (fi), L 2/(1+S) (n) c H- S (n), H 1+ "(n) c L^n) 

for < s < 1 and rj > imply that for any 9 € [1/2, 1] the force F given by (10) satisfies (6) with 
some 5 < 1/2. In the case 9 € (0, 1/2) we rely on the inequality 

Wf-gWs^CWflUvWgW^, f,g€H\Sl), 
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which holds for any < s < 1 and < a < 1 — s (see [11, Lemma 1.4.1]) and implies that 



i 

ii/i---/iii.<cnn/iiii-*. o<5<i, 
i=i 

where < 5 < 1 depends on s and I. Using the latter inequality we can find that 

iidiv{|vuirv«i-iv« 2 rvu2} 11-20 < ciiiv«irvui-ivu 2 rv« 2 iii-2e 

< C , ||ui-« 2 ||2-tf(||«l||2 l +||«2|2 l ), 

where m is either q or r. Thus (6) holds for every < 9 < 1. We also have 

n(u) = J ${u{x))dx + J \Vu{x)\ q+2 dx- J | Vu(x)\ r+2 dx — J u(x)p(x)dx, 

where <&(s) = <p(£)d£ is the antiderivative of ip. It follows from (11) that there exist 7 < \\ 
and C > such that $(s) > -7« 2 /2 - C for all s£R. This implies (7). 

In the case of the von Karman model, we have that the Airy stress function v(u) defined in 
(12) satisfies the inequality 

||[iii,i;(iii)] - [u 2 ,v(u 2 )]\\- v < C(||ni|| 2 + ||« 2 || 2 )||ui - u 2 \\ 2 - v 

for every 77 € [0, 1] (see Corollary 1.4.5 in [11]). Thus, (6) holds for every < 9 < 1 with 5 = 9. 
The potential energy II has the form 



n W = \f M u )\ 2 ~ 2([«,*b] " 2p)u] dx 



and possesses the properties listed in Assumption 1.1(F), see, e.g., [11, Chapter 4]. 

One can also see that the Berger model satisfies Assumption 1.1(F) for every < 9 < 1, and 
(6) holds with 5 = 9; for some details see [6, Chapter 4] and [10, Chapter 7]. 

Damping terms: Now we consider possible forms of the damping operator. 
Case 9 = 1: Our main example is (9) under the following conditions: 

• cr (si), <7i(si, s 2 , S3) and g(si, s 2 ) are locally Lipschitz functions of Sj 6 R, i = 1, 2, 3; 

• 0"o(si) > 0, (7i (si, s 2 , s 3 ) > and #(si,s 2 )s 2 > for all s«; 

• there exists q > such that 

- ^(Dt < c R ie - r la + I6r + \&\ q + IC2 r + I8n 

for all £ = (6,6,6), 8 = (£,8,8) e ™ch that |6|, [81 < 

• there exist (?i < 4 and q 2 < 2 such that 

i^co - 5(di < c R [(1 + i6i 91 + m qi m + (i6P + i8m& - e 2 \] m 

for all £ = (6,6), 8 = (8,8) G ^ such that 161, 181 < 
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In particular, we can take 

<ti(u, Vu) = a 10 (u) + a u (u)\Vu\ r and g(u,u t ) = go(u)u t + gi{u)u\, 

where <7io, <7n, go and gi are nonnegative locally Lipschitz functions, r > 1. 

Case # = 1/2: We consider the case when er = and o\ is independent of Vu, i.e., we 
consider a damping of the form 

D(u,ut) = -div [ai(u)Vu t ] + g(u,u t ). (43) 
One can see that Assumption 1.1(D) holds with 8 = 1/2 under the following conditions: 

• o"i(si) and 5(51,52) are locally Lipschitz functions of S{ G R, i = 1,2; 

• ci(si) > and <?(si,S2)s2 > for all 

• the function g satisfies (42) with some q\ < 3 and (72 < 2. 

We can also consider more general (anisotropic) damping operators D(u, ut) which are defined 
variationally by the formula 

(D(u,ut),ip) = / <Hjkl(u)u Xi x j t'tl>xix j dx + y2 / bij(u)u Xlt ip Xj dx + \ c(u,Vu,u t )ipdx 

for every ip € (H 2 n Hq)(Q.) under appropriate positivity and smoothness hypotheses concerning 
the coefficients. 

5.2 Wave equation with strong damping 

As an example, we can also consider the following wave equation on a bounded domain in R 3 
with a nonlocal damping coefficient: 

utt - <TQ{\\u\\ v )/±u t + <Ti{u)u t - Au + (p(u) = f(x), u\ m = 0. (44) 

We assume that rj < 1 and the following conditions concerning the damping functions gq and o\ 
are valid: (i) (Tj(s) are locally Lipschitz functions of s G R, j = 0, 1; (ii) <Jo(s) > 0, cri(s) > for 
all sGi; (iii) there exists q\ < 3 such that 

- ^(f )i < c (i + ler 1 + ir n \z - c i, e,r e r- 

We also assume that the source term ip G C 1 (R) possesses the properties 

liminf {(^(s)s" 1 } > -Ai, \<p'(s)\ < C(l + \s\ q ), s G R, <? < 4, 

|s|— >oo 

where Ai is the first eigenvalue of the Laplacian with the Dirichlet boundary conditions. 

In this case we can apply Theorem 3.1 for the model in (44) with 8 = 1. We also note that 
basing on a requirement like (8) we can also cover the case when r\ = 1 in (44). As for the case of 
the critical growth exponents (q± = 3 and q = 4) of the damping coefficient <j\ and the force ip, we 
cannot apply here our abstract approach. This case requires a separate consideration involving a 
specific structure of the model. 

In a similar way, we can also consider the wave model (44) in arbitrary dimension d and with 
another structure of the damping operator. 
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